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Abstract 

A compact general integral formula is derived from which the fermionic contribution to the 
onedoop coefficient in the perturbative expansion of the MS coupling in powers of the bare lattice 
coupling can be extracted. It is seen to reproduce the known results for unimproved naive, staggered 
and Wilson fermions, and has advantageous features which facilitate the evaluation in the case of 
improved lattice fermion formulations. This is illustrated in the case of Wilson clover fermions, and 
an expression in terms of known lattice integrals is obtained in this case which gives the coefficient 
to much greater numerical accuracy than in the previous literature. 
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I. INTRODUCTION 



When transforming results from lattice simulations into a continuum scheme such as MS 
it is often desirable to know the perturbative expansion of the renormalized coupling in 
powers of the bare lattice coupling. This is useful as an intermediate step for relating the 



MS coupling to the coupling defined in nonperturbative lattice schemes such as the ones 



based on the static quark potential 



and Schrodinger functional 



4 , and is also needed 



to translate bare lattice quark masses into the MS scheme (see, e 



The one loop 



coefficient in the expansion is of further interest because it determines the ratio of the lattice 
and MS A-parameters . [lol. [ill. [12I. [13I. \\\ . Moreover, the one loop coefficient is also 



needed for determining the two loop relation between the couplings, from which the third 
term in the lattice beta-function (governing the approach to the continuum limit) can be 



determined 



13. 



15 



16 



17|,ll8|,ll9|. 



In this paper we derive, for general lattice fermion formulation, a compact general integral 
formula from which the fermionic contribution to the one-loop coefficient in the perturbative 



expansion of the MS coupling in powers of the bare lattice coupling can be extracted. The 
motivations for pursuing this are as follows. First, given the plethora of lattice fermion 
actions currently in use, and the likelyhood of new ones or improved versions of present ones 
being developed in the future, it is desirable where possible to have general formulae from 
which quantities of interest can be calculated without having to do the calculation from 
scratch each time. Second, it is desirable to have independent ways to check the computer 
programs used these days to perform lattice perturbation theory calculations via symbolic 
manipulations. Third, by reducing the calculation to a managable number of one loop 
lattice integrals one can more easily achieve greater numerical precision than with symbolic 
computer programs. This is important, since, as emphasized in 20j, the one loop results 
need to be determined with very high precision to achieve reasonable precision in the two 
loop result. As a demonstration that the general formulae of this paper are useful in this 
regard, we apply them to obtain the fermionic contribution to the one loop coefficient in the 
case of Wilson clover fermions 2jJ to almost twice as many significant decimal places as in 
the previous literature. 

As reviewed in sect. II, determining the fermionic contribution to the one loop coefficient 
reduces to determining a constant cj arising in a logarithmically divergent one fermion loop 
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n = 2 



n = 3 



n = 4 




FIG. 1: The lattice Feynman diagrams for the gluonic n-point function with internal fermion loop 
for n = 2,3,4. 

lattice Feynman integral I (am), which has the general structure 

I {am) = ^ log(a 2 m 2 ) + c 7 (1) 

Here a is the lattice spacing and m an infrared regulator fermion mass. The numerical 
factor in the log term is universal, whereas Cj depends on the details of the lattice fermion 
formulation. I (am) arises from the one fermion loop contribution to the gluonic 2-point 
function, and it is from this that it was evaluated in previous works for specific lattice 
fermion formulations. However, Ward Identities allow I (am) to also be evaluated from the 
gluonic 3- or 4-point functions. In this paper we evaluate I (am) from the one fermion loop 
contribution to the gluonic 4-point function. In this case there are five lattice Feynman dia- 
grams to consider rather than the two diagrams for the gluonic 2-point function - see Fig. [TJ 
Nevertheless, evaluation of I (am) from the 4-point function turns out to be advantageous. 
The diagrams are evaluated at vanishing external momenta without the need to first take 
momentum derivatives, and we find three nice properties: 

(i) Only one of the five diagrams is logarithmically divergent - it is the first n = 4 diagram 
in Fig. [TJ The other four diagrams are all convergent. 

(ii) The logarithmically divergent diagram is not affected by changes in how the link vari- 
ables are coupled to the fermions (e.g., it is unchanged by adding staples, clover term etc.). 
Consequently, it is the same for improved and unimproved versions of the lattice fermion 
formulation (provided the free field formulations are the same). 
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(iii) The four convergent diagrams, or subsets of them, vanish when the lattice Dirac opera- 
tor is sufficiently simple. In particular, they all vanish for unimproved Wilson and staggered 
fermions, also when the Naik term |22|] is included. 

Thus for improved versions of Wilson and staggered fermions the only new quantities to 
compute relative to the unimproved case are the four convergent one- loop lattice integrals. 1 
The main result in this paper is a general integral formula for I (am) obtained by eval- 
uating the contributions from the five n = 4 Feynman diagrams in Fig. [1] for general lattice 
fermion formulation, from which the desired constant c/ can be extracted. Specifically, we 
do the following: (a) evaluate the contribution from the logarithmically divergent diagram, 
deriving a quite explicit general formula which is seen to reproduce previous results for the 
cases of unimproved Wilson and naive/ staggered fermions, and (b) derive formulae for, and 
describe a straightforward procedure for evaluating, the contributions from the four con- 
vergent diagrams. We illustrate this in the case of Wilson clover fermions. The general 
formulae lead to integrals to which the method of Ref. [23| can be applied, reducing the 
integrals to basic lattice integrals that are already known to high precision. The application 



of our result to other lattice fermion formulations such as Asqtad staggered fermions 



and overlap fermions 



24 ] 



251 ] will be made in future work. 



The paper is organized as follows. Sect. II reviews the one loop expansion of the MS 
coupling in the bare lattice coupling, using the background field approach. In sect. Ill we 
derive an initial expression for I (am) as the sum of contributions from the five n = 4 diagrams 
in Fig. 1 and point out the properties (i), (ii) and (iii) mentioned above. Rather than 
evaluating the diagrams directly, we infer them from perturbative expansion of the fermion 
determinant, which is easier. From this the general formulae and applications mentioned 
in (a) and (b) above are derived in sect. IV and V, respectively. The concluding sect. VI 
describes applications to be carried out in future work, as well as the possibility of deriving 
similar results for the gauge-ghost contribution to the one loop coefficient. Also, it is pointed 



out that the present results are relevant for a previous proposa 
field action on the lattice from the lattice fermion determinant 



for constructing the gauge 



261 ] . Specifically, our results 



give an expression for the coefficient of the Yang-Mills action that arises in that proposal. 



1 This is only true in the staggered fermion case when the Naik term is not used. More on this in the 
concluding section. 
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Some technical details of our calculations are provided in two appendices. 



II. GENERALITIES OF THE ONE-LOOP RELATION BETWEN BARE LAT- 



TICE COUPLING AND MS COUPLING 

The gauge field quantum effective action T(A) in Euclidean spacetime can be expressed 
(prior to gauge fixing) as 

e -r(A) = f VQV ^ e -j,S YM (A +9 Q)~ftD(A +g Q)^ (2) 



where Sym is the Yang- Mills action and D(A + gQ) is the Dirac operator coupled to A + gQ 
where A is the background gauge field and Q the quantum fluctuation field. The gauge 
group is taken to be SU(N) with the fermion fields being in the fundamental representation. 
T(A) has the loop expansion 

T(A) = - 2 S YM {A) + Y X {A) + g 2 T 2 (A) + ...+ g 2 ^T n (A) + ... (3) 
9 

Expanding T(A) in powers of A, the term quadratic in A has the form 

T W(A)= fT^(p)^A^ P )A b u (- P ) (4) 

where A a {p) denotes the Fourier transform of A a {x). Gauge invariance of T(A) and rotation 
symmetry imply that the 2-point function has the form 

r (2) G< = -5 ab ( P \ v -p, Pu )\- 2 - Mp) - 9 2{n - 1} vn(p) -■■■} (5) 



The relation between the MS renormalized coupling g and bare lattice coupling op at 



one loop can be obtained by requiring equality of T(A)j^g and T(A)i at up to one loop 



13: 



- 2 S YM {A) + T^ms = \s YM (A) + T^Ajua (6) 
9 9a 

Expanding each side in powers of A, the equality betwen the quadratic terms gives, in light 
of©, 2 

\ ~ "i(pW = 3 " Mp)iat (7) 
9 9o 



2 In the lattice theory there is only hypercubic rotation symmetry. However, this suffices to obtain ([5]) in 
the lattice setting up to terms which vanish for a — > 0. 
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The one-loop contribution to the 2-point function is given by a sum of 3 terms: gauge 
field loop, ghost loop and fermion loop. Consequently, for Nf flavors of massless fermions, 
u i(.P)ws an d vi{p)i at have the general forms 

Mp)m = -ft MP 2 /V) + cy + N fC ^ (8) 

vi (p)iat = -A) log(a Y) + cf at + N f c{ at (9) 



where fx is the mass scale in the MS scheme, a is the lattice spacing, and 

h=^(*%-*£) ( io ) 



16tt 2 V 3 J 3< 

The constant c g arises from the gauge and ghost loop contributions; it depends on N and the 
gauge fixing parameter, while & arises from the fermion loop contribution. In the lattice 
case c 9 and & also depend on whatever parameters are present in the lattice gauge and 
fermion actions. (E.g., for lattice Wilson fermions c{ at depends on the Wilson parameter.) 
From (CD and (ED (ED we get 

1 1 



-a = -2+A>log(aV)+fo (11) 
9 9o 



where 



l = l 9 + N f ll (12) 

(Our notations I ,Iq, Iq follow the papers of Panagopoulos and collaborators, e.g., 14]) The 
relation bewteen g and go up to one loop now follows from ffTT| : 

g 2 = g 2 (l - g 2 [{3 log(aV ) + 1 ] + 0(g$ ) (14) 

Al 



so, from 
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ITTT) or (|T4"D the relation between the lattice and MS A-parameters is obtained 
111: 



^ lat _ gio/2/3o 

The focus of our attention in this paper is the fermionic contribution to Iq > i- e -> m 

1« 



for general lattice fermion formulation. The continuum constant ct-^ is well-known (see, 



e.g., [16]): 



J = °_ 

^Ms 7271-2 
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(16) 



Therefore, to determine l Q we need to determine the lattice constant cf at for general lattice 
fermion formulation. 

To study c{ at it suffices to consider the Nf = 1 case which we restrict to henceforth. We 
proceed by expanding the lattice 2-point function r|f](p)^ in powers of (the components 
of) the external momentum p. Although r^(p)"^, itself is infrared finite, the terms in the 
expansion are individually infrared divergent. To deal with this we introduce a fermion mass 
m as regulator; it renders the expansion terms finite. Since gauge invariance is maintained 
when m is introduced, the expansion of the one fermion loop contribution to up to 

2nd order in p results in an expression of the general form 



J 



ns(p) 



fJtU 



5 a V<V - V^uMam) + R(p 2 /m 2 )} 



(17) 



/, l—loop 

up to terms which vanish for lattice spacing a — > 0. Here I{am) is a logarithmically divergent 
one loop lattice Feynman integral given by 



I(am) 



f, l—loop 



p=0 



(// v , no sum) 



r (2) (V) n ' 

1 lat \P) yuv 

ij-). By the structural result of |27j] it has the general form 3 

1 



(18) 



I(am) 



24vr 2 



log(a m ) + cj 



(19) 



where the constant cj depends on the details of the lattice fermion formulation. R(p 2 /m 2 ) 
in (TUD denotes the continuum limit of the remainder term after expanding [rg^J]/,!-^ 
to 2nd order in p. It is convergent by power-counting and therefore coincides 28( with the 
corresponding (known) continuum term: 



R{p 2 /m 2 



Aix 2 



24tt 2 
1 

247T 2 



p 



dxx(l — x) log x(l — x) — - + 1 



77T 



log(p 2 /m 2 



1 



47T 2 JO 

\og(p 2 / m 2 ) + R{m 2 /p 2 ) 



dxx(l — x) lo£ 



x(l — x) + 



m 



P 



Substituting (JTSJ) and O in flTTp gives 



(20) 



^laliP) 



fj,v 



S ab (p 2 5^ - p^Pu) log(a 2 p 2 ) + c/ + R(m 2 /p 



f, l—loop 



,2 /^2^ 



247T 



(21) 



3 The factor 1/247T 2 in the log term is fixed by universality 27[; it is minus the fermionic term in /3q with 
N f = l. 
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Comparing this to the expression (jSJ) for T\J t (p) ab u with V\{p)i a t given by (jUJ), we see that 

cL = Cl + R(0) (22) 

In fact -R(O) is precisely the MS constant cL^: explicit evaluation of the integral in fl20l) at 
m = gives R(0) = -5/72tt 2 . It follows from (03) that 

= -cj (23) 

Thus the issue is to determine the constant cj appearing in the logarithmically divergent 
one-loop lattice integral I (am) in (|18p - (|19p . In this paper we are going to derive a compact 
general formula for I (am) for general lattice fermions which can then be used to determine 

Cl . 



III. AN INITIAL FORMULA FOR I(am) 

In this section we derive an initial general formula for I (am). It is the starting point for 
deriving more explicit formulae in the subsequent sections. 

Let D be a general lattice Dirac operator which is translation-invariant and transforms 
covariantly under gauge transformations and rotations of the 4-dimensional Euclidean hy- 
percubic lattice. The ansatz for the link variables is 

U^(x) = e aA ^+hV (24) 

(/i=unit vector in the positive //-direction). The Fourier-transformed field A^(p) is defined 
via 

A,(x + \afi) = [ i M (p) e *-(*+l«A) . (25) 

Here and in the following J p (- ■ ■) = Jl^ a -^yr (•••)■ The continuum gauge fields are A^(x) = 
A^(a;)T a where {T a } are generators for SU(N) normalized such that tr(T a T b ) = — \5 ab . (We 
take the T a 's to be anti-hermitian, absorbing into them the imaginary unit that multiplies 
A in other notations.) 

Expanding the link variables in powers of A gives an expansion of the lattice Dirac 
operator, 

oo 

D=Y,D n (26) 

n=0 
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Translation-invariance implies that each D n can be expressed in momentum basis in the 
form 

D n (k',k) = a n ~ l S(pi + ...+p n + k-k')d n (ak\api,...,ap n )^ 1 ... fln A IJil (p 1 )---A l j in (p n ) 

Jpi,...,p n 

(27) 



The continuum limit requirements on D are 



— do{ak) i^ v ky-\-m (28) 



a 



d x (ak | ap)^ ^ 7^ (29) 

The mass m enters via an additive term am in d^(ak) (we suppress the m-dependence in 
the notation). 

To derive a general formula for I (am) it is useful to exploit the fact, clear from ([2]), that 
the one fermion loop contribution to T(A) is given by the fermion determinant: 

r(A) 7 , = - log detl) (30) 

where, in the present lattice setting, D is coupled to the link variables (j21j) of the background 
gauge field A. Decomposing D as 

D = D + D int (31) 

where Dq is the free field operator and Di nt the interaction part, we start from D = Dq(1 + 
Do l Di n t) to obtain the expansion 



log det D - log det D = Tr log(l + D l D 



int 1 



00 (— IV 

= Z] — J — TrKD^Dint) 1 ] (32) 
1=1 1 

From this, using D int = J2 n >iD n with D„ given by (l2"7j) . an expansion in powers of A is 
obtained. It has the general form 

log det D — log det D = r det(Pi> • • • ,Pn-i) m -M«* r Aui(Pi) ' • ' A ^ n (p n ) (33) 

(Pn = -(pi + • • • + p n -i))- The "n-point function" T^ t (p 1} . . . ,p n -i) m -Mn can be Taylor- 
expanded in powers of (the components of) pi,...,p n _i. Gauge invariance and mass- 
dimension considerations imply that the terms in the expansion can be combined to take 
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the form 4 

logdet D - logdet D = -I(am)S YM (A) + £ —S r (A) (34) 

r>l 

up to terms which vanish for a — > 0, where each S r (A) is a gauge invariant function with 
mass dimension r > 1. The fact that the coefficient of Sym(A) is —I (am) is inferred from 
(HU) and (J30D. 

It is clear from (|34p and (13"U1) that I (am) can be evaluated from either the 2-, 3- or 4-point 
function in (1331) . For the 4-point function the relation is 

I(am) = -Y^lt (0,0,0) V v , no sum (35) 

for any choice of \i and v with fi ^ v. Usually I(am) is evaluated from the 2-point function 
via ( jl~8]) . but we will see in the following that it is advantageous to use ( 135]) instead. To 
evaluate it we start by collecting the terms containing 4 powers of A in the the expansion 
of log det D in (|32|) . These are 

tfTrKD^D^] + tfsTrKD^D^D^Dz] 
+ i^L(2Tr(D 1 D l D 1 D,) + Tr^D, 1 D 2 ) 2 ]) + Tr(D^D A ) (36) 

Inserting the expression ( 1271) for each D n , and evaluating the traces in momentum basis, 
the T ( ^ t (p 1 , p 2 , Ps) ^ lf i 2 ^3^4. m fl32l is readily worked out. The resulting expression for I(am) 
from (1351) is as follows. From the functions d n in (I2T1) we define 

d n (k)» v ..» n = d n (k | 0, ... , 0) w ... Mn (37) 

Then 

J(am) = 1(1,1,1,1) (am) + ^(1,1,2) (ara) + ^2,2)(H + /(i, 3 )(om) + J (0 , 4 )(am) (38) 

where 

1 / >7r <i 4 &; 

^(i,i,i,i) (am) = - J J2^y tr d (k)~ 1 di(k) ll d Q (k)~ 1 dx(k) v d (k)~ 1 d 1 (k) fl d (k)~ 1 d 1 (k) v 



In order for the link variables ([24]) to transform in the correct way under lattice gauge transformations 
the transformation of the continuum gauge field A needs to be modified in the lattice setting 1C 



However, one still finds that the mass-dimension zero function in (|34|) must be ~ Sym(A) - this can be 



inferred from the argument in T. Reisz's proof of renormalizability of Lattice QCD [29(; see also §5 of [lfj 
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(39) 

j-K ^k 

^(i,i,2)(H = - / j7r-utrdo(k)~ 1 d 1 (k) ll d (ky 1 d 1 (k) l/ d (ky 1 d 2 (k) llu (40) 

J — 7T l^TTj 

1 /' 7r d 4 fc 

1(2,2) (am) = -J j^—^trd (k)~ 1 d 2 (k)^d (k)~ 1 d 2 (k)^ (41) 

rTT d A k 

I(i,3){ am ) = J j^^ trd o(k) d 1 (k) ll d (k) d z (k) vlJ , v (42) 

PTT d^k 

1(0,4) (am) = - y j—^trdoik^d^k)^^ (43) 

for any choice of // and z/ with fi ^ u (no sum over repeated indices). The traces are 
over spinor indices alone. Each integral corresponds to an n = 4 Feynman diagram in 
Fig. [Tjin sect. I and the notation reflects the structure of the corresponding diagram. The 
above integrals could have been derived directly from the Feynman diagrams, with fermion 
propagator and vertices read off from (1271) . but it is easier (and equivalent) to use the fermion 
determinant expansion as we have done above. 

The expression of I (am) as the sum of the integrals fl39l) - P3l) is our initial general formula. 
It has a number of advantageous features which we note in the following. 

In light of ( 128]) it is clear that the integral I(i,i,i,i) (am) diverges logarithmically for am — > 
while the other integrals (140]) - (143]) are all finite in this limit. The divergent integral 



necessarily 27| has the general form 



^(1,1,1,1) ( am ) = 7^"^ log(a 2 m 2 ) + a (44) 

up to terms which vanish for a — ■> 0. Then cj is given by 

cj = cj + J (1)1>2) (0) + J (2)2) (0) + J (1)3 ) (0) + 7(0,4) (0) (45) 

Next we point out that I(i,i,i,i) (am), and hence Sf, depend only on the free field lattice Dirac 
operator. This is because di{k)^ in (139]) is determined by do{k) : for any lattice fermion 
formulation we have 

d 1 {k) lt =\-£ r d {k). (46) 

This is derived in Appendix A. Therefore, a change of gauging of D (e.g., by adding staples, 
clover term etc.) affects only the finite integrals ()40l)-(l43l). 

A further advantage becomes apparent when considering the description of D in terms of 



lattice paths 



30 



31| (we will discuss the path description more explicitly in sect. V). It is 
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easy to see that the vertex function d n (k \pi, . . . ,p n )tn-vn receives contributions only from 
lattice paths which contain a lattice link parallel to the /^i-axis, preceded at some point by 
a link parallel to the /i2-axis, preceded at some point by a link parallel to the /^3-axis, and 
so on. In particular, if the lattice paths describing D are straight lines, as is the case for the 
Wilson-Dirac, naive and staggered operators, then d n (k)^...^ n vanishes unless the /ij's are 
all the same. It follows that the finite terms (l4"0]) - (l4*3]) all vanish in this case (since fi ^ v 
there). Thus, for such lattice Dirac operators, I (am) is given entirely by 1 1) (coti) . If 
staples "/i^//" are attached to such an operator then (j4T)l) - (l4"2"I) are non-vanishing while ( )43|) 
still vanishes. If a clover term (~ closed paths around plaquettes) is added then all the 
finite integrals (JHUHdSD are potentially non-vanishing. However, as we will see in sect. V, 
d^(k) v ^ u and d^k)^^ both vanish in the clover case, so the terms (03) (S3D also vanish 
there. 

These properties make (|38|) - (|45|) useful for evaluating /(am) and c/ in practice for specific 
lattice fermion formulations, especially for improved formulations. E.g., when the improved 
D differs from the original one by a more complicated choice of gauging the only new 
quantities that need to be evaluated are the finite integrals ( l4T)l) -(l4"3 l) with m = in do(k). 
The quantities d 2 (k) fMU , d 3 (k) u ^ v and d^(k) ^ Ufll/ appearing in their integrands are generally 
easy to determine in practice as we will see in sect. V. Besides these, the integrands only 
involve the free fermion do(k) and its derivative (recall f|4"6]) ). 

Taking the initial formula for I (am) in this section as starting point, we go on to derive 
more explicit expressions in the following two sections. 

IV. FORMULAE FOR (am) 

For simplicity we use the notations d = d^(k) and <9 M = g|- in the following. We assume 
that 



is a scalar, i.e., trivial in spinor space. Note that this is a free field statement; it holds for 
all lattice Dirac operators of current interest (naive, staggered, Wilson, overlap,...). 
Recalling (|4"6T) : di(k)^ = —id^do , and using the relations 



= dld — dodjj 



(47) 





(48) 
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evaluation of (|39|) leads to 

W>(»0 = -J^ — ^gl (49) 

for any choice of p, v with p ^ u , where 

Xp, = -6dj9 M (9„do ^4 ^0) + dl(d u A 4 d v d ) - <9j(A 9^4 <^d ) 

+ ^(^A 4 «9^o) - 2<9JA 9„4 0„d o (50) 

The details of the calculation are given in appendix B. To evaluate this expression further, 
we now assume that do (the free field momentum representation of the lattice Dirac operator 
D) has the general form 

d = Pa + A (51) 

where p a (k) and X(k) (which includes the mass term) are real scalar functions, and the 
gamma- matrices are hermitian, so that d = —i-f a Pa + A and Ao = p 2 + A 2 . This is the 
typical free field form of lattice Dirac operators of interest in practice; it covers the naive, 
staggered, Wilson, and overlap operators and their improved versions. For simplicity we 
make the further assumption that 

d^dydo = for p ^ v . (52) 

This holds for Wilson, naive and staggered fermions, also when the Naik term is included, 
but does not hold for the overlap operator. (The extention of the following to the general 
case d^dydo 7^ is straightforward but tedious, and we defer it to a future article where the 
specific results for the overlap operator will be derived.) 

Substituting (jSTl) into (J15]) - (J3D|) . a calculation using fl52l and the other mentioned prop- 
erties gives 

J (1)lil)1) (am) = T J ^ {p{k y + X{KamW 0* ^ ") (53) 

where 

V = dlpldlpl-i2(d,p,)\d v p v ) 2 + 2dl\ 2 dlpl-2A(d,\) 2 (d uPv ) 2 

- I2dl\ d v X p v d uPu - 4(<9 M A) 2 A d 2 u \ + 4A 2 9jA 8 2 u \ . (54) 
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From this, using integration by parts (in a similar way to the calculations in appendix B), 
a more compact expresion can be obtained: 

Wam) = - £ 0- t ( «W + _ _ i{BA log Ao)2 ) (55) 



(recall A = p 2 + A 2 ). Although this expression is more compact, the preceding one f l53|) -fl54 
seems more useful for evaluating /(i,i,i,i)(am) in practice. 

In the remainder of this section we show that these general formulae readily reproduce 
the previously known results for the specific cases of (unimproved) naive, staggered and 
Wilson fermions. (Recall from sect. Ill that I (am) is given entirely by I(i,i,i,i)(a?7T.) in these 
cases.) 

A. Naive and staggered fermions 

For a naive fermion, p a = sin k a , A = am and fl54l) reduces to 

Yfj,v(k) = 4 cos 2k ll cos 2k v — 12 cos 2 k^ cos 2 k u (56) 

up to terms which are 0(a) for a — > 0. Since a naive fermion decomposes into 4 degenerate 
staggered fermions [32.], I (am) for a staggered fermion is obtained by dividing the naive 
fermion result by 4. The invariance of (l56l) under k a — *■ k a + 7r allows the integration domain 
in (153p to be restricted to [— tt/2, 7r/2] 4 at the expense of an overall factor 16. Thus we find 

W 2 d^k o cos 2A; U cos 2A;^ — cos 2 /c u cos 2 k u , , 

/ am = JV t / ——7 ^ ^ 57 

J-^(2nr ((am)2 + E CT sin 2 A; a ) 

up to terms which vanish for a — > 0, where N t is the number of fermion "tastes" (16 for a 
naive fermion, 4 for a staggered fermion). This is precisely the expression for I (am) derived 



previously in Eq. (6.8) of 12| where the contribution to the gluonic 2-point function from a 
naive/ staggered fermion loop was evaluated. (I(am) corresponds to 11^ in [l^.) 

After changing variables by k a — > k a /2 the integral (f5T|) can be evaluated by the method 



of Ref. 



231 ] . It expresses the integral in terms of certain basic lattice integrals that were 



evaluated numerically to high precision in [23]. We find the following result: 

I (am) = N t [-^—(\og(2am) 2 + lE - F ) + ^-Z ] (58) 
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where 7# is the Euler constant and Fq , Zq are numerical constants listed to high numerical 



accuracy in Table 1 of 



231 ] . Thus the divergent part has the correct universal structure, and 



the constant cj (= — Iq ) that we are after is given by 

-cj/N t = 1 (jr - 7j5 -log(4))-^-Z 



24tt 2 V" u " ov "V 144" 

0.0026247621012431485... 



(59) 



In Ref. 



12J this constant was denoted P4 and our value agrees with the one in Eq. (6.16) of 



that paper. We have obtained it to much higher numerical precision here though, thanks to 
the high accuracy to which F and Z were evaluated in [231 ] . 



B. Wilson fermions 

In this case p a = sin k a , A = r X)<*(1 ~ cos k a ) + am where r is the Wilson parameter, and 
5^1) reduces to 



Y^ y (k) = 4 cos 2k^ cos 2k v — 12 cos 2 k^ cos 2 

+ r 2 ^8 sin 2 fc^(cos 2/c„ — 3 cos 2 fc„) — 6 sin 2k /1 sin k^ cos A;^ + 16 cos 2k ii cos ^ sin 2 

a 

+ r 4 ( — 8 sin 2 fc^ cos fc„ sin 2 (A; Q /2) + 16 cos k^ cos A;^ sin 2 (A; a /2) 



up to terms which are 0(a) for a — > 0. Substituting this into (153]) we recover the previous 
result of Kawai et al. for the one fermion loop contribution in Eq.(3.24)-(3.25) of Ref. 10]. 
(Dimensional regularization of the infrared divergence was used there, but the result is easily 
transformed into mass-regularized form (cf. §5.1 of 23j) and then coincides with ours.) The 
constant Cj is related to the constant L defined there by 



(61) 



The integral I (am) can again be evaluated by the method of Ref. 23J; in fact this was done 



in that paper for r = 1, and the expression for L in terms of the basic constants is given in 
Eq.(6.8) of 23]. From that, /£ is obtained to high precision in the Wilson fermion case: 



it = -cj 



0.00669599933173308... 



(62) 
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The Wilson fermion case was also considered independently by P. Weisz in ll|. Note 
that for Wilson fermions (d^p^) 2 + (c^A) 2 = sin 2 + cos 2 k^ = 1, hence fl55l) simplifies to 

Ham) = - jT -0^ ( ^ - \(d,d v log A ) 2 ) (63) 

thus reproducing Weisz's expression in Eq. (13) of [11]. The constant l[, denoted as P3 there, 
was evaluated numerically but to low precision - a systematic error of ~ 2% was estimated, 
and this is indeed the case when comparing the value for P3 in Eq.(22) of H[ with the high 
precision value in (1621) above. 



V. EVALUATION OF THE FINITE INTEGRALS 

In this section we derive general formulae for the finite integrals (j4"0]) - ()43p and describe 
how these can be straightforwardly evaluated, using staple and clover terms as illustrations. 
Evaluation of I (1,1,2) (am) . Starting from (|40p . calculations of the same type as in appendix 
B lead to 

^ {am) = -Lwr a 2 (64) 

(for any choice of /i, v with p 7^ v) with the notations d = do(k) , A = dl(k)do(k) and 
d 2 fj, v = d 2 (k) fiu . Specializing as before to the case where d has the general form d = i^ a p a +X 
and satisfies d fl d u d = for p, ^ v, and assuming that the gamma-matrix structure of d 2iW 
has the general form 

d 2 (k)^ = j^e^k) (p ^ v) (65) 

(El leads to 

J (1)1>2) (am) = ~J n j^y 4 (p^TW ( } 

(There are no terms involving derivatives of e tiu (k) since these all cancel out.) 
Evaluation of I( 2 . 2 )(am) . From ( I4TT) we obtain 

T , , 1 f n d 4 k tr d\d 2llv d\d 2iiV 

l m (am) = -J_^ M (67) 



Specializing as before, this leads to 



/■Tr dfLfc \ 2 e 2 

I( 2 , 2) (am) = -2J_^ (2?r)4 ( ^ + ^ 2)2 (68) 
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Evaluation of In^(am). From (jl2l we find 

J (li8) (am) = y_ )r ^ ? ^ (69) 

with the notation = dz(k) vlxl ,. Specializing as before, and assuming that dzvpv has the 
general form 

dz{k) vlu/ = 7 M e^ - z/^ (70) 



691) leads to 



J (1 , 3) (am) =4^^ (71) 

Evaluation of (am) . From ( H3i) we find 

J ^ (am) = ~L(2^ A (72) 
with the notation d^ vllv = d^k)^^. A specialized formula can be worked out as in the 
previous cases (we omit the details). Note that terms in d^ u ^ v involving a product of two 
or more gamma matrices give vanishing contribution to the trace in (!72|) . 

To evaluate the integrals in practice one needs to determine dz^k)^ , d^(k) vlxv and 
d4$) nvya, (fi 7^ v) for the lattice Dirac operator D. This can be done straightforwardly 
from the description of D in terms of lattice paths, as we now describe. In the path descrip- 
tion, D is expressed as 

Dip(x) — - cpTpUp[x, x + aA-p] ip(x + aA-p) (73) 
a j, 

where the sum is over a collection of translation-equivalence classes V of oriented lattice 
paths. Each has an associated numerical constant c-p and element Tp of the Clifford algebra 
generated by the gamma-matrices. Associated with each equivalence class V is a vector 
Ap with integer components: it is the difference in lattice units betwen the start- and end- 
points. Up[x, x + aA-p] denotes the product of the link variables along the representative 
path for V starting at x + aAp and ending at x. 

The Clifford algebra- valued functions d n (k\pi, . . . ,Pn)ti v -nn m the expansion (j2"6"I) - ([27]) of 
D are found by expanding the link variable products Up [x, x + aA-p] in (1731 in powers of 
the continuum gauge field A. In the present case we only need to determine these functions 
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at vanishing "external momenta", i.e., c/ n (A;) Ml ... Aln = d n (k\0, . . . , 0) A11 ... Atn . Therefore we can 
take the continuum gauge fields {A^} to be constants. Then Up[x, x + aAp] is independent 
of x, and its expansion in powers of the continuum gauge field has the form 

n 

The functions dnik)^...^ are then given by 

^ra(&)/ii-'Mn = dn(k)m-i*n (75) 

where 

(*)/*■■ = c ^<-/,„ eiA "' fc (76) 

Thus, given the description (!73l) of D in terms of lattice paths, the problem of determining 
d n (k) n-L-fMn is reduced to determining the coefficients t^...^ in the expansion (jTIj) of the 
link variable product Up with constant A, which is generally straightforward in practice. In 
the present case things are further simplified since we only need to know d n {k) ^...^ when 
fix ■ ■ ■ fi n is \w , v^iv , and Since \x ^ v we can replace each link variable U^ 1 appearing 

in the product U-p by 

Uf 1 ^ { (77) 
[ 1 if a <£ {//, i/} 

To illustrate this, consider the case of a 'V/zz/ staple" f/p = UyU^U' 1 . The relevant parts 
of the expansion are found by 

U V U,U- 1 - (l + 4,)(l + A M )(l-4,) 

= - AvApA,, + other (78) 

where "other" refers to terms which are not proportional to A fl A u , A u A fJi A u or A^A U A^A U 
and therefore do not contribute to d2(k) liv , d^{k) nu ^ v or d^(k) . Thus the relevant con- 
tributions from the staple are 

d^(k)^ = -cpT v e lk - , ^(fc)^ = - C7 ,r P e^ , <(^W„ = (79) 

where CpTp are whatever factors that accompany the staple Up in the expression (173|) for 
the lattice Dirac operator. 
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As another example we now consider the Sheikholeslami-Wohlert clover term for Wilson 



clover fermions 



2l|: 



C — C sw — ^ GapPa, 



(80) 



a/3 



where c sw is a tunable constant, a a p = §[7^, j v ], and P a p{x) is a sum of products of link 
variables around oriented plaquettes in the a/3-plane starting and ending at x. The relevant 
contributions in this case come from the [iv parts: 



(81) 



For constant gauge fields, Pn U is explicitly given by 



P 



u^u- l u,u„ 



u-'u^u- 1 



(82) 



Expanding this as described above, the relevant expansion coefficients are found to be 



1 j Uj/^i/ 



0, and u^vpv = 0. Since the plaquette paths are closed we have A-p = in all 



cases. Using this, it follows from f!75l) - fl76l) that 



<k{k)^ = -c sw \-i i a v , d 3 (k) Vflu = , d A {k) 







(83) 



independent of k. Note that, as discussed in sect. Ill, there are no contributions from the 
Wilson-Dirac operator since in its path description the paths are all straight lines. 

The finite integral contributions to I (am) can now be determined in the case of Wilson 



clover fermions. By (1831) and the previous formulae, /(i,3) 



'(0,4) 



0. The other integrals 



are determined by substituting = sinfc M , A = rE ff (l — cosk a ) + am and e Mi/ = —c sw /2 
into fl66l) and (1681 . Ta king the Wilson parameter to be r = 1 and evaluating the integrals 
by the method of Ref. 23J we find, in the a — > limit, 

2 



1(1,1,2) (0)/c s 



and 



W0)/cr 



.2 



3 -^(l,0)-^(2,-l) 
0.005046714024535753066... 

-^(2,-2) 

-0.02984346719542684815... 



(84) 



(85) 
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where .F(l, 0) , J- (2, — 1) and J- (2, —2) are certain basic convergent one loop lattice integrals 



231 ]. whose numerical values are given to high precision in Table 2 of 



defined in §4 of Ref. 
that paper. 

Recalling (|45|) . collecting the results for Wilson clover fermions with r = 1 we have 



(( 1621 means the numerical constant given in Eq. ( 1621) . etc.) This agrees with the previous 
literature but gives the numerical constants to much greater precision. The previously most 
precise values were those in Eq. (14) of Ref. [l?| where (1841) and (I85I) were given to 11 decimal 
places. We have obtained them here to 20 decimal places, thanks to the high precision with 



see {4] and the references therein. 

VI. CONCLUDING REMARKS 

The focus in this paper has been on deriving the general integral formulae for I(am), 
confirming their correctness by checking that they reproduce the known results in the cases 
of staggered, Wilson and clover fermions, and in the process developing general techniques 
for evaluating the formulae. In doing this we were able to express I(am) in those cases in 
terms of basic one loop lattice integrals that have already been evaluated to high precision 
in 23j. This had already been done in 23) in the Wilson fermion case, but the results for 
the staggered and clover cases (( |59|) and (|84l-( l85i) . respectively) are presented here for the 
first time. The further applications of the general formulae are left for future work, and in 
the following we discuss some possibilities for this. 

Obvious targets for future applications are the various improved versions of staggered 
fermions. These formulations involve "smearing" of the link variables to reduce flavor sym- 



which the basic integrals in Ref. 23] were evaluated. For earlier results for these quantities, 



metry violations; specifically there is the "Fat- 7" link |33|] and HYP link [34]]. Since these 
differ from unimproved staggered fermions only by a choice of gauging, evaluating I(am) 
should be straighforward: expand the relevant products of link variables in powers of the 
constant continuum gauge field to determine d 2 (k) fiu , d 3 (k) u ^ u and d^k)^^ as described in 
sect. V; then, from the formulae in sect. V the conver gent integrals I(i,i,2), A 2 : 2 ); ^(1,3) an d 



-^(0,4) can be explicitly evaluated by the method of Ref. 23j ; this will express the integrals in 
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terms of basic lattice integrals that have already been calculated to high precision. Adding 
these to the already known /(y^Jaro) (= I (am) for unimproved staggered fermions) then 
gives I (am) in these improved cases. 

Of more interest, however, is the case of 0(a 2 ) improved "Asqtad" staggered fermions 24] 
that are currently used by the MILC collaboration to generate the ensembles used in high 
precision lattice simulations 35]. Besides smearing of link variables, this formulation also 



contains the Naik term 



221 ] which modifies the free field staggered Dirac operator. Therefore, 



1,1,1) i s n °t the same as for unimproved staggered fermions in this case, and the method 
and results of Ref. 23j do not apply. I(am) may still be readily determined in this case from 
the general formulae and techniques of this paper, but it will be necessary to numerically 
evaluate the one loop lattice integrals that arise. (The approach of Ref. [36] could be used 



for this.) The relation between the MS coupling and bare lattice coupling has already been 
determined to 2 loops for Asqtad staggered fermions via symbolic computer program 
Determining the fermionic contribution to the one loop coefficient independently via the 
(semi-)analytic approach of the present paper would provide a useful check on the computer 
program. 

Another future application is to overlap fermions. The fermionic contribution to the one 



loop relation between the MS coupling and bare lattice coupling in this case was calculated 



via symbolic computer program in 



in 
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131 ]. and the two loop relation was subsequently calculated 



19]. Application of the results of the present paper will allow the one fermion loop 



contribution to be obtained from numerical evaluation of one loop lattice integrals without 



3j_in this way will be 



191 ] also relies on. It 



the need for symbolic computer programs. Reproducing the result of_ 
a useful check on the computer program, which the two loop result 
will also allow the one fermion loop contribution to be calculated more easily and to higher 
precision for any value of the overlap parameter that one wishes to consider. I emphasize 
that, in the approach of the present paper, the technical problem of calculating the one 
fermion loop contribution for overlap fermions is greatly reduced compared to the usual 
approach followed in 13J. The simplification comes about because the determination of the 
quantities d 2 (k) IJ , v , d 3 (k) u ^ u and d^(k) ^ Ufll/ in the expansion of the overlap Dirac operator 
can be done with constant continuum gauge fields. Expansion of the overlap Dirac operator 
in powers of constant continuum gauge fields is relatively straightforward, and has already 



been successfully used [37] to reproduce results obtained via symbolic computer program in 
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another context 381. 5 



Having found useful general formulae for the fermionic contribution to the coefficient 



in the one loop relation between the MS and bare lattice couplings, a natural question is 
whether a similar approach is possible for the contribution from the gauge and ghost loops. 
In fact this seems quite possible. In the background field approach the gauge field in the 
action is A + gQ, and it is the part quadratic in the quantum fluctuation fields Q that 
determines the one loop contribution to the effective action. Coming from the functional 
integral of a quadratic term, it is clear that this can be expressed as a functional determinant, 
both in the continuum and lattice settings. There is also the Faddeev-Popov determinant 
through which the ghosts arise. One can then expand the determinants in powers of A as 
done for the fermion determinant in this paper. The constant to be determined in this case, 
Ciat m ©; can again be found from the quartic term in the expansion (via Ward Indentities). 
In the fermionic case considered in the present paper, the possibility to introduce an infrared 
regulator mass was crucial for deriving the general formulae. This can also be done in the 
gluonic-ghost case: note that Q transforms under gauge transformations by Q — > UQU^ 1 
so a regulator mass term mtrQ 2 may be introduced without breaking gauge invariance. 
This is currently being pursued, and I hope to be able to present general formulae for the 
gluonic-ghost contribution to the one loop coefficient in future work. The hope is that this 
may allow allow a (semi-) analytic calculation of the one loop contribution in the case of 
improved lattice gauge actions, which has recently been evaluated via symbolic computer 
program in 

Finally, the results of this paper are relevant for a previous proposal for constructing the 
gauge field action on the lattice from the lattice fermion determinant [26|]. Set 77 = am 
and regard m = rj/a as a function of 77 and a. Expanding r^(pi, . . . ,i?n-i)^i-/u« i n f|33|) in 
powers of the momenta, without taking the a — > limit, leads on dimensional grounds to 
the following variant of ([34]) : 

logdetDfo) - logdet D ( V ) = -I{rj)S YM (A) + £> r <? r (A,7/) (87) 



Ref. 



371 ] is the work mentioned in the concluding section of 38j . 
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where S r (A, rj), a function of A and 77, has mass-dimension r. Now taking a — > we obtain 6 

lim ( logdet D(ri) -\ogdet D (rj)) = -I(ri)S YM (A) (88) 

(Strictly speaking this depends on the sum in ( 1H71) being convergent, which requires the 
continuum gauge field A to be sufficiently weak.) Thus the coefficient of the YM gauge 
action obtained from the lattice fermion determinant is seen to be which can be 

evaluated from the general formulae in this paper after replacing am by 77. 7 In the rj — > 
limit we have 

m «° -j^a M^) + c ' ( 89 ) 

Although calculation of the constant cj was the main focus in this paper, the general formulae 
can just as well be used to determine I(rj) for general values of 77. The integrals will need to 
be evaluated numerically for the specific values of rj that one considers though. 
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APPENDIX A: DERIVATION OF THE FORMULA (0BD FOR di(fc)^ 

We derive this from the path description (1731 -(|76j) of the lattice Dirac operator. From 
fP76|) we have 

a%{k) = c v T v e lAv - k (Al) 

= opTrtfe**** (A2) 

6 Since to = 77/a, taking a — * with 77 fixed amounts to taking a simultaneous continuum and large mass 
limit. The fact that the large mass limit of the lattice fermion determinant effectively gives a lattice gauge 
action was discussed before in [3^]. 

7 The relationship between I{r)) and the coefficient c s (r)) in [26| is I(rj) = 2c s {rj). 
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Thus to show the relation (jlBjl : di(k)^ = —id^doik) it suffices to show 

< = (Ap) m (A3) 

From the definition (1741) it is easy to see that counts the number of links of the path V 
that lie along the /^-direction, counted with sign ± depending on whether they are oriented 
in the positive or negative /i-direction. But this is precisely the /z-coordinate of the difference 
(in lattice units) between the start- and end points of (a representative path for) V, i.e., the 
/^-coordinate of A-p. Thus we have found that flA3j) holds. 

APPENDIX B: DERIVATION OF THE GENERAL FORMULA FOR 1(1,1,1,1) (am) 

The general formula (l49|) - (!50l) is obtained by evaluating the integrand in the initial ex- 
pression ( 1391) as follows. For notational simplicity we omit the "tr", write d, A for do, Aq, 
and use "~" to denote equality up to terms which vanish upon taking the trace or terms 
which are total derivatives and therefore give vanishing contribution to the integral. 

d^d^d dT x d v d dT^d^d d~ x d v d 
= d^d^dydd^d^dyd 

= i^dt _ M^^ay _ 

111 

- ~£?, 9 ^ 9 v d 9 ^ 9 » d + 9 n (^2 ) d] 9 » d 9 ^ 9vd ~ ^ ( ) ®^ ^ dyd ^ ® ud 

~ ^ 9 ^ d 9 ^ 9 » d ) + ^3 ( 9 l A ^ 9 - d ^ 9 » d + 9 » A 9 ^ d] ^ S 9 » d )) 

(Bl) 

The second term here is re-expressed using 

$ d u d S d u d = (d u A - d v S d)S d v d = d u A d ] d v d - A d v S d v d (B2) 

to get 

\d u (^)dlAJd u d-^dlAd u Sd u d ~ -^ 2 (^d v {dlASd v d)-dlAdJd v d) 

(B3) 

The third term in ( 1B1I) is re-expressed as 

^{^d^S d v dd) d v d) ~ -^cfrtfdvdSdvd) (B4) 
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The resulting expression for the integrand, given by the first term in (1B1[) plus (1B3|) plus 
Ml . leads to the claimed formula (H9D-(|50D. 
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